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Electron collisions with O2 at scattering energies below 1 eV are studied in the fixed-nuclei ap- 
proximation for a range of internuclear separations using the ab initio molecular R-matrix method. 
The ^Ilg scattering eigenphases and quantum defects are calculated. The parameters of the reso- 
nance and the energy of the bound negative ion are then extracted. Different models of the target 
that employ molecular orbitals calculated for the neutral target are compared with models based on 
anionic orbitals. A model using a basis of anionic molecular orbitals yields physically correct results 
in good agreement with experiment. An alternative method of calculation of the R-matrix is tested, 
where instead of performing a single complete diagonalization of the Hamiltonian matrix in the 
inner region, the system of linear equations is solved individually for every scattering energy. This 
approach is designed to handle problems where diagonalization of an extremely large Hamiltonian 
is numerically too demanding. 
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I. INTRODUCTION 

After nitrogen, molecular oxygen is the most abun- 
dant molecule in the Earth's atmosphere. Therefore, its 
study is crucial to our understanding of planetary at- 
mospheres, while also providing useful insight into the 
physics of gaseous discharges and laboratory and astro- 
physical plasmas. Theoretical and experimental research 
into electron interactions with O2 has attracted signifi- 
cant scientific attention. Since a complete summary of 
the relevant references on this topic exceeds the scope of 
the present study (see the recent review by Itikawa [ij 
and references therein), we mention only those directly 
related to the problem addressed here. An experimental 
study of the resonant vibrational excitation of O2 by the 
low-energy electron impact by Linder and Schmidt 
shows oscillatory structures in the cross sections due to 
the ^Ilg resonance. Celotta et al. Q later measured the 
electron affinity of O2 using molecular photodetachment 
spectroscopy. Those two complementary studies 0, Q 
provide a deeper understanding of the structure and dy- 
namics of the ^Ilg state of 

Several previously published theoretical treatments 
based on the J?-matrix method deal with the ^Ilg res- 
onance in the fixed-nuclei approximation . All those 
calculations yield the energy of the ^Ilg resonance at the 
equilibrium geometry of O2 more than 0.6 eV above the 
value obtained from the experimental spectra 0, H, 
These studies are unable to reproduce the bound state of 
at larger internuclear separations without further ad- 
justments of the fi- matrix poles [1, • Higgins et al. [8] 
used the adjusted results of the ah initio calculations to 
study resonant vibrational excitation by electron impact 
and compared the computed cross sections with the ex- 
perimental work of Field et al. fll| . The resonant vibra- 
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tional excitation of O2 at low energies has been revisited 
by [13 ■ Their theoretical treatment of the nuclear motion 
is based on the "boomerang model" and that utilized the 
same adjusted positions and widths of the resonances in 
the fixed-nuclei approximation calculated by Noble et al. 

The energies of the ^Ilg resonance and the bound state 
of were later calculated more accurately by Ervin 
et al. [1^ using the methods of quantum chemistry de- 
signed originally for bound states. However, to our 
knowledge, no reliable ah initio calculation exists that 
provides reliable scattering phase shifts at low energies 
(< 2eV) for different symmetries. 

Significant scientific attention has been paid to elec- 
tron collisions with O2 at energies above 2eV as well. 
Teillet-Billy et al. found in their study based on an 
extended version of multichannel effective range theory 
that the ^n. resonance dominates the electronic excita- 
tion of the ^Ag and ^11+ states. This process was later 
studied experimentally [l4l - [l6j and the results confirmed 
the existence of the resonances predicted by the previous 
calculations 

Good correspondence between theoretical and exper- 
imental results for the excitation of the ^Ag and ^Ej" 
electronic states by electron impact has been achieved 
(see paper and references therein). However, the re- 
maining discrepancies between the theoretical 0, fisj 
and experimental cross sections [3[l7, 18] for excitation 
of the Herzberg pseudocontinuum (^S^ , '^Au, "^S^) still 
constitute a challenge 

More recently, Tashiro et al. [l^ presented another 
_R-matrix calculation of the elastic and electronically in- 
elastic electron collisions with O2 at energies above 5 eV 
that shows better agreement with the experimental re- 
sults than previous ab initio studies. 

The main goal of the research presented here is to de- 
velop an advanced fi-matrix calculation of electron colli- 
sions with O2 in the fixed-nuclei approximation at scat- 
tering energies below 1 eV. In particular, improvements 
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over previous calculations are sought to provide a more 
physical energy and width of the ^Ilg resonance and en- 
ergy of the anionic bound state than was obtained in pre- 
vious i?-inatrix studies. The scattering eigenphases and 
the quantum defects discussed here can be employed in 
the investigation of the resonant vibrational excitation 
of O2 based on the non-local resonance model [20|or 
energy-dependent vibrational frame transformation [2l|. 
The treatment of bound and continuum states of in- 
troduced here will be later adapted to study other sym- 
metries of the anionic complex that are relevant to O2 
photodetachment and to the vibrational dynamics of that 
process. 

In comparison with previously published ab initio 
studies, a more physical representation of the anionic 
complex electronic structure studies is achieved through 
the use of a more advanced treatment of electron correla- 
tion and polarization in the inner region and by employ- 
ing a basis set of molecular orbitals optimized for the neg- 
ative ion instead of the neutral target. Note that in other 
contexts, this idea has also proven to be beneficial [2^ . 
Current implementations of ab initio electron-molecule 
scattering theory rely on expansions of the wave func- 
tions associated with both the neutral iV-electron tar- 
get and the {N + 1) -electron system which use the same 
truncated basis of molecular orbitals. Although orbitals 
optimized for the neutral target are traditionally used in 
molecular R-matrix calculations, it is far from obvious 
whether that is universally the most appropriate choice. 
The present study compares i?-matrix calculations per- 
formed using optimized neutral target molecular orbitals 
with those performed using optimized anionic orbitals, in 
order to ascertain which set is more appropriate. 

Another goal of this study is to test the feasibility of 
an alternative method of calculation of the R-matrix that 
has the potential to treat larger basis sets and more ex- 
tensive configuration interaction. The well-established 
approach preferred in the UK ff-matrix codes requires 
a single full diagonalization of the Hamiltonian matrix 
in the inner region [23l . HJ] . The R-matrix is then easily 
calculated for any scattering energy, since its dependence 
on the energy is extremely simple in the eigenrepresen- 
tation of the Hamiltonian. The development of more 
efhcient _R-matrix methods capable of treating molecules 
having a larger number of active electrons is one of the 
goals driving this research direction. The requirement of 
a full Hamiltonian diagonalization in the present imple- 
mentation has made it quite challenging for the J?-matrix 
method to handle more complex polyatomic molecules 
and molecules for which electron correlation plays a key 
role. Evaluation of all the Hamiltonian matrix elements 
and their storing in the computer memory prior to the 
full diagonalization is necessary in the present implemen- 
tation. That approach is computationally much more 
demanding than the more economical methods of the 
quantum chemistry of the bound states, where only a 
small fraction of the spectrum is calculated and the itera- 
tive methods of the diagonalization are routinely utilized. 



Those require only evaluation the matrix elements that 
are used in the actual iteration. In addition, the com- 
plete set of the eigenvectors that have a dense structure 
is necessary in the present implementation of the ab initio 
molecular f{-matrix method. That further raises issues 
with the memory limitations. 

The method presented in this study employs the so- 
lution of a system of linear equations having the same 
dimension as the Hamiltonian matrix, individually for 
every scattering energy. A number of other methods 
have previously been introduced to overcome this diffi- 
culty in solving problems requiring very large basis sets. 
The mo st p romising among them is the partitioned R- 
matrix [24 . [25j . It consists of a single calculation of 
few lowest eigenvalues and eigenvectors of the Hamilto- 
nian matrix and of a model-like approximation of the 
rest of the spectrum. While the partitioned R-matrix 
method retains the advantage of a single diagonalization 
of a (reduced-dimensionality) Hamiltonian matrix, the 
alternative method presented here requires solution of 
the linear system of equations for every scattering en- 
ergy. However, it is free of any model-like assumptions, 
and therefore, is a complementary approach to the par- 
titioned R-matrix method. It is hoped, however, that 
in combination with generalized quantum defect meth- 
ods or multichannel effective range theory, the quantities 
to be calculated will be comparatively smooth functions 
of energy and can accordingly be calculated on a coarse 
energy mesh. Note that the method of calculation of 
the R-matrix by solving the system of linear equations 
has been formulated previously by Collins and Schneider 
[2^. Their linear-algebraic approach utilizes the static 
exchange approximation of the electron interactions in 
the inner region. 

The ab initio theoretical description of electron col- 
lisions with O2 at low scattering energies is challenging. 
Any successful treatment of the complicated polarization 
effects [13, and electron correlation in the inner region 
requires a large basis set of configurations. Depending on 
the details of the model, the usual approach based on the 
full diag onalization of the inner-region Hamiltonian ma- 
trix |23l H^l can become intractable. That makes this 
system a suitable candidate to test the performance and 
limitations of the alternative method investigated here. 
As is discussed below, the linear solution method proves 
to be feasible and advantageous for calculations where 
the dimension of the Hamiltonian matrix in the inner re- 
gion exceeds 40000. Several different models of the tar- 
get are introduced below and, where possible, the perfor- 
mance of the traditional full diagonalization method is 
compared with that of the direct linear solution method. 

The rest of this paper is organized as follows: Sec. [H] 
describes the method adopted to calculate the R-matrix 
without any diagonalization, by direct solution of a linear 
inhomogeneous system of equations. Different models 
of the neutral target are discussed in Sec. IIIII and the 
construction of the {N + l)-electron Hamiltonian matrix 
is described in Sec. lIVl The results are presented in Sec.lVl 
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and the conclusions are summarized in Sec. I VII 

Atomic units are used throughout, unless stated oth- 
erwise. 



II. R-MATRIX BY LINEAR EQUATION 
SOLUTION (LES) 

The idea of the R-matrix method is the solution of 
the Schrodinger equation within a finite reaction volume 
n of the configuration space. The scattering properties 
of a many-particle system are known once the normal 
logarithmic derivative of the wave function and relevant 
surface amplitudes are specified on the surface E enclos- 
ing the reaction volume. The goal of the theory is to 
determine this information in the form of an R-matrix. 
The reaction volume in the molecular ab initio fi-matrix 
method is specified by a sphere of radius vq chosen such 
that ri < ro, Ti being the distance between the ith elec- 
tron and the center of mass of the molecule, ro is large 
enough to contain all the complicated interactions of the 
electrons within the inner region, and sufficiently large 
that the target wave functions in the open and weakly- 
closed channels are negligible in the outer region. Interac- 
tion of the scattering electron with the target in the outer 
region is well approximated by the long-range multipole 
and dispersion potentials that couple different scattering 
channels [IJ . The total wave function can be written as 


^^(r,^)^^-*,^^^,^^, (1) 

1=1 

where r is the radial coordinate of the scattering elec- 
tron and w denotes all other spin-space coordinates of 
all the electrons (including the spin and angular coordi- 
nates of the scattering one). includes the electronic 
state of the target as well as the spherical harmonic of 
the scattering electron, Fip{r) is the scaled radial wave 
function of the scattering electron in the outer region 
(r > ro) in the channel i. Nch denotes the total number 
of scattering channels retained (both open and closed), 
index /3 = 1 , . . . , Nch denotes different linearly indepen- 
dent solutions for the total energy e of interest. If the 
radial derivative of Fipir) is denoted as F-p{r), then the 
K-matrix on the surface E is defined 29] as 

R^{FF:-')r. (2) 

and is calculated by solving the Schrodinger equation in 
the inner region. Then it is used to match the solutions in 
the outer region and to calculate the iC-matrix or other 
quantities that characterize the scattering process (23l . 

The approach to calculation of the _R-matrix tested in 
this work is based on the noniterative variational formu- 
lation of the fi-matrix method introduced by Ref. [131 
that has been used extensively in a number of problems 



(see Ref. 0] and references therein). Solutions of the 
time-independent Schrodinger equation at energy e in the 
inner region obey 

H^p = e^p, (3) 

where H is the electronic Hamiltonian of the [N + 1)- 
electron system and e is the total energy. \I'^(r, oj) can 
be expressed using a set of real basis functions yfc(r, w) as 
^^(r, w) = yfc(r, a;)Cfe/3. Each of the basis functions 
j/fc can be expanded on E as 

yk{ro,uj) = — $j(a;)uife(ro). (4) 

2=1 

As is well known from Robicheaux [3l| and Greene and 
Kim [s^], the i?- matrix can be expressed as 

Rij ^^Uik{^^'^)kk>Ujk', (5) 

kk' 

where 

Tfefc' =2 / ykie-H-L)yk'dV. (6) 
Jn 

The integration is performed over the reaction volume 
fl. Notice that the matrix F is symmetric due to the 
presence of the Bloch operator L [2^ |3^ . The basis 
set used in the UK J?-matrix program suite allows for 
a close-coupling expansion of the total {N + l)-electron 
wave function in the inner region 

^'^(r,w) = i^C,fe/3i$,(w)M,fc(r) 

i.k 

+ Y.Dppx'^^'M. (7) 
p 

where Uik{r) are the radial parts of the continuum or- 
bital introduced to represent the scattering electron in 
the inner region and their values at ro are in general 
non-zero. The angular parts are included in $i(w). The 
choice of the continuum orbitals depends on the symme- 
try of the target electronic states. These two are cou- 
pled to give the correct overall spin and spatial symme- 
try of '^fj{r,ijj). Index i denotes the scattering channel 
in the outer region and characterizes the electronic state 
of the target as well as the partial wave of the scatter- 
ing electron. Therefore, each basis function that appears 
in the first sum in Eq. ([7|) has non-zero amplitude on 
the boundary E and is associated with one scattering 
channel. Furthermore, the electrons must obey the Pauli 
principle and they are anti-symmetrized by operator A. 
The second summation in Eq. ([7]) involves antisymmet- 
ric {N + l)-electron configurations Xp^^ that have zero 
amplitude on the boundary E and where all the electrons 
occupy the orbitals associated with the target (L^ con- 
figurations [1^). Notice that in the terminology of the 
variational _R-matrix method l29l-[3ll the basis functions 
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in the first sum in Eq. ([7]) correspond to the open part 
of the basis, while the second summation corresponds to 
the closed part. 

The eigenstates of the target and the {N + l)-electron 
basis functions are both expressed in terms of the com- 
plete active space configuration interaction (CAS CI) 
[34| . It is worth mentioning at this point that the dimen- 
sion of the " closed " part of the basis (second summation 
in Eq. ([7])) is typically much larger than the dimension 
of the "open" part that has non-zero amplitudes on the 
boundary E. 

The modified Hamiltonian matrix for the inner region 
= iL + L. calculated using the basis expansion ([7]) is 
evaluated in the UK R-matrix codes as well as the sur- 
face amplitudes Uik{ro). The matrix F defined by Eq. [6] 
can be easily calculated for each scattering energy e of 
the interest and Eq. (O can be used to calculate the R- 
matrix. 

The product {r~^}^j,/ Ujk'{ro) is implemented as 
a solution of the linear system of equations, which is the 
most computationally demanding step of the calculation. 

A celebrated aspect of the approach to the calcula- 
tion of the _R-matrix available in the UK _R-matrix suite 
nowadays is that the matrix _ff q is diagonalized only once 
and the ener gy d ependence of the ff-matrix is calculated 
analytically [231 0, [3l|. However, the complete set of 
eigenvalues and eigenvectors in the given basis set is nec- 
essary for accurate evaluation of the _R-niatrix using the 
expansion in the eigenstates. Beyond certain size of the 
basis set the matrix storage hits the memory limit or the 
time necessary to diagonalize the modified Hamiltonian 
becomes too long for practical calculations. 

The approach based on solution of the linear system for 
each individual energy of the interest becomes favorable 
in those cases, since both the time and memory required 
to solve one system of linear equations is significantly 
smaller than the time required to completely diagonal- 
ize a matrix of the same size. Furthermore, while exist- 
ing computer routines for complete diagonalization are 
usually based on the full matrix storage, several mod- 
ern computer implementations of state of the art lin- 
ear solvers are based on more economical sparse stor- 
age schemes. Since Hr, is often a sparse matrix, the 
approach introduced above enables calculations that use 
larger CI models than the diagonalization-based method. 
A more efficient parallel implementation of the linear 
solvers than of the algorithms for the complete diago- 
nalization makes the method presented above even more 
favorable for large-scale calculations performed on high- 
performance computer clusters. A calculation of the R- 
matrix for a single value of the energy using Eq. ([5]) can 
be executed in the parallel mode in those environments 
and calculations for multiple energies can be trivially par- 
allelized. 

Several other alternatives to the complete diagonaliza- 
tion of have been proposed and used [11-25, 35, 36]. 
They are based on the accurate calculation of a modest 
number of the lowest eigenvalues and eigenvectors, while 



the rest of the spectrum is approximated using various 
models. These approximate approaches preserve the ad- 
vantage of the single diagonalization and of the analytical 
energy-dependence of the _R-matrix. The most promising 
method among those is the partitioned _R- matrix ,25i] . 
The method described above is not based on the diago- 
nalization of and is free of any model-like assump- 
tions. 



III. MODELS OF THE NEUTRAL TARGET 

All the J?-matrix calculations presented here were per- 
formed using the polyatomic UKRmol program suite 
[23I HJl , which uses a basis set of Gaussian- type orbitals 
(GTOs) in the inner region. The irreducible represen- 
tations of the D2h point group are used, as this is the 
largest abelian subgroup of the true Dock symmetry point 
group of O2. The notation of the Dooh point group is 
used throughout the rest of this paper, unless otherwise 
stated. 

The results are presented for a range of internuclear 
separations from i?„ = 1.9 a. u. to i?„ = 3.5 a. u., in- 
side of which the potential energy curves for the ground 
electronic states of both O2 ('^^g ) and (^Ilg) reach 
their minima. The fixed-nuclei scattering phase shifts 
and smooth quantum defects obtained for that range are 
essential for future theoretical studies of resonant vibra- 
tional excitation by low-energy electrons and these en- 
tities also arise in the theoretical description of vibronic 
coupling in photodetachment. The molecular orbitals 
optimized using the state-averaged complete active space 
self-consistent field (SA-CASSCF) method implemented 
in the program package MOLPRO [13, HI] are employed 
in all the fi-matrix calculations presented here. 

The neutral target is represented by one of two dif- 
ferent sets of the molecular orbitals. The first is calcu- 
lated using an extensive GTO basis set of atomic natural 
orbitals (ANO) ^39,] and the SA-CASSCF molecular or- 
bitals are optimized for the neutral molecule (three lowest 
electronic states are averaged with equal weights). This 
selection of the primary GTO basis set is motivated by 
Ref. [2^, where it was successfully used to calculate the 
electron affinity of 02- Furthermore, Jones and Tennyson 
[40| found that this GTO basis is necessary to reproduce 
within CAS CI models the static dipole polarizabilities of 
various diatomic molecules containing oxygen, although 
the case of O2 is not discussed in that study. Previously 
published fi-matrix studies [3, [3 suggest that polariza- 
tion effects play an important role in the electron colli- 
sions with O2. The ANO basis set is expected to repre- 
sent these effects more accurately than the more compact 
Gaussian basis sets employed by the previously published 
_R-matrix calculations. Two different models of neutral 
O2 based on the ANO GTO basis and on the neutral 
molecular orbitals were tested. Two additional models 
of the target are introduced below, where Dunning's cc- 
pVTZ [41] GTO basis set is employed. The convergence 
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of the scattering calculations with respect to the size of 
the CAS and the quality of representation of the correla- 
tion and polarization effects can be estimated from com- 
parison of the scattering eigenphases and the energies of 
the stable negative ion calculated using these different 
models. 

The dominant valence electronic configuration of the 
ground state of O2 is 

It is natural to include the orbitals 3f7„, 27rg and 2tTu in 
all the CAS models, since their occupation numbers in 
the ground state are larger than 0.01. All of the models 
presented here consist of 8 active valence electrons. The 
CAS models considered in the previously published R- 
matrix studies d, 0, SI] include all 12 valence electrons 
and a smaller set of the active orbitals than the CAS 
models introduced here. However, our preliminary tests 
suggest that the correlation effects due to the electrons 
2(7 g and 2au can be neglected for the range of the collision 
energies considered here. 

The first CAS model can be expressed as 

(lo-g2o-glcr„2o-„)^(3crg3cr„l7r„27r„37r„l7rg27rg)®. (8) 

The _R-matrix calculations based on this CAS include 
the two energetically lowest target states from each ir- 
reducible representation that does not contribute to the 
static dipole polarizability of the ground state. Tashiro 
et al. suggested that the polarization effects were not 
sufficiently represented in previously published i?-matrix 
studies d, 0| ■ As a first step towards the improvement of 
this deficiency, the 30 lowest states from the irreducible 
representation ■^Hu and the 42 lowest states from the irre- 
ducible representation (both have a non-zero dipole 
coupling with the ground state) are included in the ex- 
pansion of the total wave function in the inner region. 
This set of target states and the CAS expressed by Eq. (|S]) 
is denoted as Model 1 in the text below. 

In order to evaluate the convergence of the scatter- 
ing calculations with respect to the number of the active 
molecular orbitals, another model of the target is intro- 
duced. The orbital 37r„ included in Model 1 is replaced 
by the orbital Aag. This CAS can be expressed as 

(lfTg2fT3lfT„2cr„)^(3o-g4(Tg3(T„l7r„27r„l7rg27rg)^. (9) 

The R-matrix calculations based on this model include 
the 40 energetically lowest target states in every symme- 
try (both singlet and triplet spin states). This model is 
denoted as Model 2 in the text below. 

The number of molecular orbitals from every irre- 
ducible representation included in the treatment of the 
inner region for every CAS model introduced here is sum- 
marized in Table m 

The scattering eigenphases calculated using the target 
models introduced above provide insight into the role of 
the excited electronic states and higher molecular orbitals 



TABLE I. The number of molecular orbitals in every irre- 
ducible representation included in different CAS models of 
the target. The last line shows the number of continuum or- 
bitals (COs) that is the same for all the R-matrix calculations 
discussed here. The correspondence between the irreducible 
representations of the point groups -D2/1 and Dooh is shown as 
well. 



Symmetry (I>2h) 


ag 


h2u, bsn 


big 


biu 


bag, b2g 




Symmetry {Dooh) 


Og, 5g 










5u 


Model 1 


3 


3 





3 


2 





Models 2,3 


4 


2 





3 


2 





Model 4 


4 


3 





3 


2 





COs 


37 


21 


18 


21 


18 
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in electron collisions with O2. However, they do not yield 
physically correct energy of the ^Hg resonance and fail to 
describe the bound state of O^. In order to solve this de- 
ficiency two additional target models (Models 3 and 4) 
based on different MOs are introduced. They employ 
the Dunning cc-pVTZ GTO basis 41] and the CASSCF 
molecular orbitals that are optimized for the ^Hg ground 
electronic state of O2 ■ The choice of the cc-pVTZ basis is 
motivated by the experience with the K-matrix calcula- 
tions based on Models 1 and 2 discussed below and by the 
intention to use the obtained results in the prospective 
calculations of the nuclear dynamics using the vibrational 
frame transformation method [21I, . The jR- matrix cal- 
culations based on the ANO GTO basis require a rather 
large J?-matrix sphere (rg = 16 a. u.), which causes nu- 
merical difficulties with the vibrational frame transfor- 
mation method. Models 3 and 4 yield spatially more 
compact molecular orbitals that can be confined inside 
the sphere with the radius tq = 10 a.u. As is discussed 
below, the fi-matrix calculations based on Models 3 and 
4 yield at the scattering energies below 1 eV more physi- 
cal results than the calculations utilizing Models 1 and 2. 
It should be kept in mind that the molecular orbitals cal- 
culated for do not have any straightforward physical 
interpretation for the internuclear distances, where the 
anion is not bound. They are used as the basis functions 
in the CI expansion of the bound target states and in the 
expansion of the {N + l)-electron scattering wave func- 
tion, where the correct boundary condition is guaranteed 
by the continuum orbitals and by the the Bloch operator 
defined on the boundary of the reaction volume. 

Model 3 consists of the same configurations as Model 2 
(see Eq. ^) and 33 target states from every irreducible 
representation are used in the expansion of the wave func- 
tion in the inner region. The most complex CAS model 
constructed in the present study (Model 4) allows the ac- 
tive electrons to occupy both orbitals Acjg and i-Ku along 
with 4cr„. It can be expressed as 

(lcrg2crglcr„2cr„)^(3crg4crg3cr„4cr„l7r„27rK37rul7rg27rg)^ 

(10) 

and it is introduced to study the stability of the _R-matrix 
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calculations based on the anionic molecular orbitals with 
respect to the extension of the CAS. The expansion of the 
wave function in the inner region includes the 33 lowest 
target states from every irreducible representation (sin- 
glet and triplet spin states) in the _R-matrix calculation 
using this model. The dimension of H^^ for this model 
is large and its full diagonalization, implemented as a 
standard method of calculation of the _R-matrix in the 
UK codes [2^, [2J would be numerically very demand- 
ing. This suggests that Model 4 is a suitable candidate 
for demonstrating the alternative approach discussed in 
Sec. In] which employs no diagonalization. 



TABLE II. Energy of the ^EJ ground electronic state of O2 
(a.u.) and vertical excitation energies (eV) for the lowest 
few excited states. The excitation energies calculated using 
Models 1-4 are compared with the experimental values quoted 
in the reference [9[. All the values are for the internuclear 
separation Re{02) = 2.3 a.u. 



ANO 



cc-pVTZ 



State 


Model 1 


Model 2 


Model 3 


Model 4 


Ref. [9] 




-149.8557 


-149.8516 


-149.8186 


-149.8430 




^A, 


1.014 


1.008 


1.005 


1.012 


0.98 




1.871 


1.828 


1.815 


1.822 


1.65 




5.875 


5.834 


5.844 


5.851 


6.12 


^A, 


6.141 


6.094 


6.099 


6.124 


6.27 




6.303 


6.254 


6.252 


6.282 


6.47 




9.458 


9.503 


9.477 


9.322 


9.25 


^A„ 


11.869 


11.926 


11.920 


11.755 


11.8 



The energy of the "^E" ground state of O2 calculated 
for the equilibrium internuclear separation Re{02) = 
2.3 a.u. is for Models 1-4 compared in Table HIl Since 
in Models 1 and 2 the neutral target is represented us- 
ing a larger ANO GTO basis and the molecular orbitals 
are optimized for the ground electronic state of the neu- 
tral molecule, it is not surprising that these models yield 
lower energy of the ground state than Models 3 and 4. On 
the other hand, the highest energy of the O2 ground state 
obtained from Model 3 is a consequence of the smaller 
GTO basis set and of the fact that the wave function of 
the neutral target is expanded in the truncated basis set 
of the orbitals optimized for O2 ■ 

The vertical excitation energies for the lowest eight 
electronic states calculated for Models 1-4 are also com- 
pared in Table HIl In general, they are in good agreement 
with each other. Note that the excitation energy of the 
lowest excited state ^Ag is for all the Models 1-4 close to 
the experimental value quoted in the reference j^. The 
agreement with the experiment is less convincing for the 
higher excited states, although the correspondence be- 
tween Models 1-4 is obvious. Table HIl as well as Table II 
in the reference [i^l suggests that the change of the pri- 
mary GTO basis as well as the presence or absence of 
the orbitals 4ag and 37r„ in the CAS do not dramatically 
affect the vertical excitation energies. 



Since one of the goals of the present study is to pro- 
vide the fixed-nuclei scattering eigenphases and energies 
of the anionic bound states for a future study of the vibra- 
tional dynamics, it is important to represent the target 
correctly also for different geometries than the equilib- 
rium one. The vibrational energy aje(02) is a suitable 
quantity that suggests how well the different models in- 
troduced above represent the potential energy curve of 
the ground state near the equilibrium. In fact all the 
Models 1-4 yield values similar to each other and to the 
previously published experimental results [Tsj . Compar- 
ison of Model 1 and Model 3 with the experimental value 
is shown in Table IIIII 

As is discussed below, the polarization effects play im- 
portant role in the scattering at low impact energies. It 
is their representation that requires rather large num- 
ber of the target states included in the expansion of the 
scattering wave function in the inner region. Thus, it is 
interesting to calculate the static dipole polarizability of 
the ground state of the neutral molecule to estimate, how 
well different models introduced above represent these ef- 
fects. The components of the tensor of the static dipole 
polarizability a of the ground state tpo can be calculated 
using the sum-over-states formula 



k>0 



r,s e {x,y,z}, 



where ipk are the excited states, dr^s are the Cartesian 
components of the operators of the transition dipole mo- 
ments and AEk is the excitation energy from tpg to (fk- 
The summation should be performed over the complete 
set of the eigenstates (only the states of the symmetry 
■^IIu or "^S" have a non-zero contribution) including the 
continuum states. Here it is performed over all the target 
states included in the expansion of the scattering wave 
function in the inner region. The only non-vanishing 
components of a for the homonuclear diatomic molecule 
are the diagonal ones azz and Uxx = cxyy One can see 
in Table IIIII that Models 1 and 3 yield a similar value 
of azz and it is more than 68% of the value that was 
previously accurately calculated by Hettema et al. [27| . 
On the other hand, both models yield significantly un- 
derestimated value of axx that is more than one order of 
magnitude lower than the accurate calculations of Het- 
tema et al. In fact, all the Models 1-4 introduced 
above yield very similar values of a and they show rather 
poor representation of the component that is perpendic- 
ular to the internuclear axis. It is worth mentioning that 
our test calculations (not published here) that utilize the 
pseudocontinuum orbitals in the CAS in addition to the 
molecular orbitals [1^ H^, SI, El] yield a higher value 
of axx- However, the R- matrix calculations with those 
huge basis sets are computationally too demanding to be 
practically performed (see also Appendix IB|) . 
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TABLE III. The equilibrium internuclear distance Re, vibrational frequency cje of the neutral target and of the ^Ilg anion 
and the electron affinity Ea- Calculations based on Model 1 and Model 3 are compared with the experimental results. The 
non-zero components of the tensor of the static dipole polarizability of the neutral target {a^x, ctzz) at the equilibrium geometry 
calculated using the same models are compared with previously published theoretical values. 







O2 {x^E- 


") 










Model 1 


Model 3 


Previously published 
values [Ref.] 


Model 1 


Model 3 


Experimental 
values [Ref.] 


Re (a.u.) 


2.30 


2.32 


2.28 fl3] 


2.55 


2.58 


2.55 [3] 


uje (eV) 


0.192 


0.185 


0.196 [13] 


0.135 


0.133 


0.137 [3] 


Oxx (a.u.) 


0.08 


0.22 


7.12 [27] 






Cfzz (a.u.) 


9.31 


9.2 


13.51 [27] 








Ea (eV) 








-0.408 


0.375 


0.448 [13] 



IV. SCATTERING CALCULATIONS 

The continuum basis for the inner region is constructed 
in the polyatomic UKRmol program using additional 
GTOs with the centers that coincide with the center of 
the R-matrix sphere. A sufficient number of these GTOs 
is diffuse enough to have non-zero values on the boundary 
of the inner region. The exponents are optimized using 
the program GTOBAS [46] and the resulting functions 
are orthogonalized on the set of the molecular orbitals. 
This procedure yields a set of continuum-type orbitals in 
the inner region. All the R-matrix calculations presented 
here include the continuum orbitals with orbital angular 
momenta / — 0,1,2,3,4,5. Their number in every irre- 
ducible representation is identical for all the Models 1-4 
and is listed in Table HI 

Models 1 and 2 are based on the ANO GTO basis that 
yields quite diffuse molecular orbitals. The correspond- 
ing _R-matrix calculations require quite a large sphere 
with radius tq = 16 a.u. A similarly large fi-matrix 
sphere was also necessary in the E-matrix studies [H, \^ 
of electron collisions with other molecules having a siz- 
able polarizability. The GTO basis set cc-pVTZ used in 
Models 3 and 4 yields target orbitals that are more com- 
pact and can be confined within the sphere of the radius 
To — 10 a.u. 

The total number of target states (summed over all 
the irreducible representations and spin states) included 
in the expansion of the scattering wave function is listed 
in Table IIVI for Models 1-4 along with the dimensions 
of the corresponding Hamiltonian matrices H_q. Quite a 
large number of target states are needed to achieve con- 
vergence of the scattering ii'-matrix and its eigenphases 
for all the models discussed here. That, in combination 
with the large number of the active electrons and orbitals, 
is the reason for the substantial amount of CPU time re- 
quired for the evaluation of all the elements of Hn (see 
the column tc in Table |IV|. Note that the CPU time 
required for the complete diagonalization of Hp {to in 
Table HVl is for Model 2 and 3 comparable with the CPU 
time necessary for the construction of Hp. The large size 



TABLE IV. The number of the target states Nt included in 
the close-coupling expansion of the total wave function in the 
inner-region and the dimension of is compared for Models 
1-4. tc is the CPU time in hours required for the calcula- 
tion of all the matrix elements, to is the CPU time in hours 
necessary to diagonalize Hn using the LAPACK subroutine 
DSYEVD and is the CPU time in hours necessary to cal- 
culate the R-matrix for a single energy using the linear solver 
PARDISO. 



Model 


Nt 


Dimension of 


tc (h) 


to (h) 


tL (h) 


1 


128 


26168 


17.7 


6.4 


0.6 


2 


640 


22456 


3.25 


3.9 


0.92 


3 


528 


22512 


4.5 


4.17 


1.2 


4 


528 


127012 


97.7 




62.2 



of the Hamiltonian matrix for Model 4 dramatically com- 
plicates the complete diagonalization of the Hamiltonian 
matrix and for this case only the method based on solv- 
ing the linear system of equations was applied. The CPU 
time required for the solution of the linear system of equa- 
tions is smaller than the time necessary for the complete 
diagonalization (see the column in Table WW\ for all 
models where both methods can be compared. However, 
it must be kept in mind that the linear system has to be 
solved for every scattering energy of interest. This makes 
the formulation of the outer-region problem in terms of 
the analytical quantum defects more favorable, since a 
smooth dependence on the energy permits the fi-matrix 
calculations to be performed for a smaller number of en- 
ergy grid points followed by interpolation when possible. 
Also, it should be kept in mind that the repetitive linear 
equation solution required at different energies is trivially 
parallelizable. Although the CPU time necessary to cal- 
culate the R-matrix for Model 4 using the linear solver 
is rather high, the efficient parallelization allows for cal- 
culating the ff-matrix for a single value of the scatter- 
ing energy on a computer with 8 CPUs in less than 12 
hours. The large CPU time required for the construc- 
tion of Hp, in Model 4 suggests that this is the most 
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time-consuming step of the J?-matrix calculation inde- 
pendently of the complete diagonalization. 

The large number of {N + l)-electron wave functions 
needed to achieve convergence for all of the Models 
1-4 considered here is in fact the usual complication 
of the ab initio calculations of the electron collisions 
with molecules with large polarizability [HI. A useful 
computational method developed to treat this situation 
efficiently is the Jl-matrix with pseudostates (RMPS) 
[231 13] , where the large number of true electronic states 
of the target is replaced by a smaller set of pseudostates 
[isj . We comment, however, that the straightforward ap- 
plication of this approach to the present problem has not 
simplified the calculations and it does not improve the 
results. Further details are discussed in Appendix IB] 

The extensive number of the target states considered 
in the inner region problem can lead to very time con- 
suming i?-matrix propagation in the outer region, if all 
the scattering channels are included in the outer-region 
problem (Eq. ([T|)) as well. This situation is similar to 
the RMPS calculations, where the R-matrix propagation 
in the outer region usually requires more CPU time than 
the complete diagonalization of Hn [3, SH] . Fortunately, 
the treatment of the wave function in the outer region can 
be simplified. As one can see in Table HI the threshold 
energy of the lowest electronically excited channel (^Ag) 
is at 1 eV above the ground state of O2. That is the up- 
per limit of the energy range considered in this study. 
All the scattering channels associated with the excited 
target states are strongly closed and although they play 
an important role in the inner region, they can be safely 
neglected in the outer region. 

Since previous i?-matrix studies [1^, H^] of the polar- 
ization effects suggest that their representation in the 
inner region plays a more important role than the polar- 
ization potential in the outer region, the J?-matrix prop- 
agation in the outer region is skipped in all the calcula- 
tions presented here and the R-matrix calculated at rg is 
used to match directly to the regular and irregular radial 
wave functions of the free particle in the outer region. 
Their linear combination determines the if-matrix and 
the scattering eigenphases. The calculations including 
different numbers of partial waves show that the ^Ilg elec- 
tronic state of for the incident electron energies below 
1 eV can be sufficiently well-represented in the outer re- 
gion by the single partial wave d, i.e. I = 2. Therefore, 
the problem in the outer region problem can be reduced 
to a single scattering channel. It is worth mentioning at 
this point that neglecting the long-range interaction of 
the target with the incident electron in the outer region 
that is predominated by the polarization potential, leads 
to a modification of the threshold behavior of the phase 
shift. While the dependence of the phase shift on the 
momentum of the incident electron k is is in the present 
results ri{k) cx for fc — 0, if the potential oc r~* was 
taken into account in the outer region, this would not be 
the leading term in the effective range expansion j49|] . 

The corresponding phase shift ri{E) for the internu- 



clear separations i?„, where the electronic state (^Ilg) 
is not bound, can be parametrized by the Breit-Wigner 
formula 



r]{E) = tan 



r/2 



(12) 



where F and are the width and position of the reso- 
nance, respectively and E = fc^/2 is the kinetic energy 
of the incident electron. The background phase shift 
770 is a smooth function of the energy and it can 
be parametrized by a low-order polynomial. This Breit- 
Wigner fitting formula is only valid for resonance energies 
well-separated from threshold, i.e. by many widths F. 

One of the goals of the present study is to provide the 
fixed-nuclei data required for the future calculations of 
the resonant vibrational excitation of O2 and photode- 
tachment of O^. To this end the analytical quantum 
defects [E) are calculated as a function of the scatter- 
ing energy for a set of the internuclear separations. The 
analytical quantum defect fio{E) is a scattering quantity 
similar to the scattering eigenphase ri{E) in the sense that 
it specifies the linear combination of the general solutions 
of the Schrodinger equation in the outer region matching 
the boundary condition © . These solutions are rescaled 
to remove the Wigner threshold factors and the analyti- 
cal quantum defect is a function of the scattering energy 
that is smooth across the thresholds and it is well defined 
for both open and closed channels [10, [Ml . The smooth 
character of the analytical quantum defects even in the 
vicinity of the resonance makes this parametrization of 
the scattering wave functions particularly favorable in 
the context of the energy- dependent vibrational frame 
transformation [2l|, [4^ . 

The reduction of the problem in the outer region to 
a single scattering channel and to a single partial wave 
I = 2 (for incident electron energies below 1 eV) implies 
that the expansion ([T]) reduces to a single term. The 
scaled radial wave function of the scattering electron can 
be for r > ro parametrized as 



Fn{r)=M[fl>{k,r) cos(7r/ 



r sm 



r>ro 



(13) 



fi{k,r) = y/2pnk ''rji{kr) is the regular scaled ra- 
dial wave function of the free particle, gfir) — 
^ 2 1 Tik^^"^ rni (kr) is the irregular scaled radial wave func- 
tion, k is the momentum of the incident electron, ji and 
ni are the regular and irregular spherical Bessel functions, 
respectively and M is the normalization factor. 

Formulation of the problem in the outer region in terms 
of the smooth quantum defects ^jp is also favorable for 
the calculations of the anionic bound states in the fixed- 
nuclei approximation. For the range of nuclear geome- 
tries, where the anionic state is bound, its energy E\, can 
be calculated by solving the equation [s^] 

TT^°{Eb) = TT - arctan \{-2Eby+^ E^ < 0, (14) 



9 



3.5 
3 

2.5 
2 

1.5 
1 

0.5 






Model 1 
Model 1b - 
Model 2 
Model 3 - 
Model 4 - 



0.2 



0.4 



0.6 
E(eV) 



0.8 



1.2 



FIG. 1. (Color online) The ^Ilg phase shift as a function of 
the incident electron energy calculated for the internuclear 
separation Re{02) ~ 2.3 a. u. Results for the Models 1-4 are 
compared. 



where the single partial wave and single target electronic 
state is assumed. The advantage of this method com- 
pared to the widely used matching of the R-matrix to 
the spherical Hankel functions (see the review [23| and 
references therein) is that both sides of Eq. are usu- 
ally smooth functions of energy and the complications 
due to the poles of the R-matrix can be avoided. 



RESULTS 



The Ilg Resonance at Equilibrium Geometry 
of the Neutral Target 



Fig. [T] shows the ^Ilg phase shift for the equilibrium 
internuclear distance of the neutral target i?e(02) ~ 
2.3 a. u. calculated using the models introduced in 
Sec. mil All the curves clearly show a narrow resonance 
with relatively small background. Model 1 yields the 
highest resonance position among all {Er — 0.838 eV). 
The expansion of the wave function in the inner region 
includes the 30 lowest target states from the irreducible 
representation ^IIu and the 42 lowest target state from 
the irreducible representation since those contribute 
to the polarizability of the ground state of the target. 
Only the two lowest target states are included from each 
of the other irreducible representations that have zero 
dipole coupling with the ground state of the target. 

Further tests showed that the higher target states from 
these irreducible representations are also important for 
achieving convergence of the phase shift within the space 
of configurations The phase shift calculated includ- 
ing the 33 energetically lowest target states from ev- 
ery symmetry in the inner region is plotted in Fig. [1] 
and denoted as Model lb. Inclusion of these additional 
states decreases the resonance position by 0.081 eV to 



Er = 0.757eV and adding even more excited states does 
not considerably change this value. 

The representation of the polarization effects by the 
CAS CI model, and the question of how well it is char- 
acterized by the static dipole polarizability of the tar- 
get, has been the subject of several studies [H, [H, |45| . 
The difference between the phase shifts denoted as Model 
1 and Model lb in Fig. [1] demonstrates that the static 
dipole polarizability is not a sufficient measure of the po- 
larization effects in the electron collisions with O2 at low 
energies. The decrease of Er with an increasing num- 
ber of excited target states was also reported by Tashiro 
et al. [iS], although that model takes into account only 
a significantly smaller set of target states; achieving con- 
vergence at low energies was not the goal of that study, 
in any case. 

In Model 2, the electrons can occupy the orbital Aag 
instead of the orbital 'inu included in Model 1. Although 
the expansion of the wave function in the inner region 
includes 40 excited states from every irreducible repre- 
sentation, the phase shift shows the resonance at slightly 
higher collision energy {Er = 0.797 eV) than Model lb, 
as one can see in Fig. [TJ 

In general, the comparison of the phase shifts calcu- 
lated using Models 1 and 2 suggests that presence or ab- 
sence of the molecular orbitals Aag and 37r„ in the CAS 
does not dramatically affect the parameters of the reso- 
nance. The energy of the ^Ilg resonance calculated for 
i?e(02) using Models 1 and 2 agrees well with the previ- 
ously published i?-matrix calculations by Higgins et al. Q 
{Er = 0.754eV). The CAS used in that study is smaller 
than in Models 1 and 3 and only the eight energetically 
lowest target states are included in the expansion of the 
wave function in the inner region. Therefore, none of 
the target states with Ilg or n„ symmetry are included 
in that study (see Table HI)) . although these states have 
significant contribution to the polarization effects. The 
absence of the higher molecular orbitals and excited tar- 
get states is compensated by the virtual orbital that can 
be singly occupied only in the (A^-l-l)-electron wave func- 
tion [7| (the configurations in the close-coupling ex- 
pansion [23| or in Eq. (O). The use of virtual orbitals 
raises the issues with the unbalanced treatment of the 
correlation in the target and in the scattering complex. 
That can lead to an ambiguity in the energy of the reso- 
nance. The models introduced in Sec. Illll are free of this 
complication. 

The resonance energy Er and width F calculated for 
the internuclear separation i?e(02) are summarized in 
Table |Vl It shows that Models 1 and 2 yield a resonance 
width F having the same order of magnitude. These 
values agree well with the results previously published 
by Higgins et al. and by Noble and Burke 6]. Our 
testing J?-matrix calculations performed with even larger 
CAS than those of Models 1 and 2 (not published here) 
show that the position and width of the ^Ilg resonance 
does not considerably change, when both orbitals 37r„ 
and 4(Tg are included in the CAS. 
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TABLE V. Energy Er and width T of the ^Hg resonance in O, 
at the equihbrium internuclear distance of the neutral target 
-Re(02). Results of the molecular R-matrix calculations using 
Models 1-4 are compared with the results in the previously 
published references. 





Er (eV) 


r (eV) 


Model 1 


0.838 


0.043 


Model lb 


0.757 


0.034 


Model 2 


0.797 


0.038 


Model 3 


0.154 


3.6 X 10"^ 


Model 4 


0.169 


1.1 X 10"^ 


Higgins et al. [7] 


0.754 


0.031 


Noble and Burke [6] 


0.700 


0.026 


Derived from experiment [9] 


0.090 


8.5 X 10"^ 



Although the results obtained using Models 1 and 2 
presented in Fig. [T] and Table |V] exhibit encouraging 
agreement with the previously published _R-niatrix calcu- 
lations 0] , they do not agree well with the experimen- 
tal study by Linder and Schmidt that found the energy 
of the ^Ilg resonance to be i^i- O.leV. A similar value 
was calculated by Ervin et al. [isj using the stabilization 
method. This suggests that the relatively good agree- 
ment of the phase shifts calculated using Models 1 and 2 
points towards the very slow convergence with respect to 
the number of the molecular orbitals included in the CAS. 
The authors of the previously published ab initio studies 
0, S 113 attribute the discrepancy between the theoret- 
ical and experimental results at energies below 1 eV to 
the insufficient treatment of the polarization effects, par- 
ticularly in the outer region. Models 1 and 2 discussed 
above employ a quite extensive GTO basis set (ANO) 
that includes a subset of diffuse functions designed to 
represent the polarization effects. It is reasonable to ex- 
pect that Models 1 and 2 account for these effects more 
than in any previously published _R-niatrix study and the 
rather large i?-matrix sphere (rg = 16 a.u.) should hold 
most of these effects in the inner region. However, even 
this improved treatment is not sufficient to provide more 
physical parameters of the ^Ilg resonance. 

The results suggest that the CAS CI representation 
employing the SA-CASSCF orbitals optimized for the 
neutral O2 is not sufficient for the reliable scattering cal- 
culations in the energy range considered here. This con- 
clusion is understandable in view of the previously pub- 
lished ab initio studies of the adiabatic electron affinity 
Ea of 02- It is another essential quantity that charac- 
terizes the ^Ilg state of at the equilibrium internu- 
clear distance of bound (i?e(0^) w 2.6 a.u. [l3|). 
Gonzalez-Luque et al. [l^ compared the value of Ea cal- 
culated using different CI models with the experimental 
value and found that the methods based on the CAS ap- 
proach do not yield even a correct sign of Ea ■ According 
to that study, really extensive multi-reference configu- 
ration interaction (MRCI) calculations are required to 



obtain a value of Ea comparable with the experimental 
results. Furthermore, Stampfuf5 and Wenzel [52] studied 
the contributions of the single and double excitations of 
the reference Hartree-Fock determinants to the binding 
energy of O2 (^Ilg) and compared it with the contribu- 
tion of the triple and quadruple excitations. Their results 
show that both contributions are comparable. Therefore, 
the set of configurations, where several electrons are ex- 
cited into the lowest few molecular orbitals (included in 
the CAS models introduced here), is not sufficient to yield 
the correct value of Ea. That requires inclusion of con- 
figurations where at least one electron occupies one of the 
higher molecular orbitals with an orbital number up to 
12 in every irreducible representation. These orbitals are 
not included in the CAS models employed in this study 
and their further extension would lead to an extremely 
demanding construction of the Hamiltonian matrix Hn- 

Since the equilibrium internuclear separation of O2 is 
only 0.3 a.u. smaller than the equilibrium geometry of 
O^, the same mechanisms are responsible for the slow 
convergence of the resonance energy Er at the equilib- 
rium geometry of 02- 

In order to improve the results provided by Models 1 
and 2, the CAS Model 3 has been introduced. It employs 
the SA-CAS MCSCF molecular orbitals of (^Hg), as 
is described above. The corresponding phase shift cal- 
culated for the internuclear separation i?e(02) = 2.3 a.u. 
is also plotted in Fig. [1] It clearly shows a sharp res- 
onance with a relatively smooth background. Fitting 
to the Breit-Wigner formula ([H]) yields Er = 0.154eV. 
This value agrees with the experimental results Q much 
better than Models 1 and 2 (see Table W} . The calcula- 
tion in the inner region includes 33 lowest target states 
from every irreducible representation and further increase 
of that number does not considerably change the phase 
shift. Interpretation of the improvement due to replace- 
ment of the orbitals optimized for the neutral target by 
the anionic molecular orbitals is not straightforward. In 
general, the orbitals optimized for the anion have more 
diffuse character than those calculated for the neutral 
target. Therefore, it is reasonable to expect that they are 
more suitable to represent the polarization effects than 
the molecular orbitals of the neutral target. This can par- 
tially compensate the absence of the higher orbitals in the 
CAS models mentioned above. On the other hand, the 
orbitals calculated for the neutral target are more suit- 
able to represent the ground and excited states of the 
target than the anionic orbitals used in Model 3. It is 
possible that the lower energy of the resonance calculated 
using Model 3 is not only a consequence of the improved 
treatment of the interaction between the target and the 
scattering electron, but to some extent also an artifact of 
less accurate model of the target. In other words, part 
of the reason of the lower resonance energy in Model 3 
is the increase of the target ground state energy with 
respect to Models 1 and 2 (see Table HT]). The CAS CI 
expansions of the target states and of the (A^-l-l)-electron 
wave function in the same truncated set of the molecu- 
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lar orbitals show different convergence with respect to 
the number of orbitals and this convergence depends on 
their character. This general complication of the ab initio 
i?-matrix calculations does not have a universal solution 
and the particular choice of the molecular orbitals appar- 
ently cannot yet be automated, but rather needs to be 
physically motivated. 

The CAS Model 3 yields the vertical excitation ener- 
gies of the target and the energy of the resonance that 
are in reasonable agreement with the experimental val- 
ues (see Table In] and Table |V| . It should be emphasized 
that these results are not adjusted by an artificial over- 
correlation of the (A^-|- l)-electron system by introducing 
virtual orbitals [1, 0, H, Ell • 

The CAS Model 4 allows the electrons to occupy both 
4(Tg and StTu molecular orbitals that are included sep- 
arately in different CAS models introduced previously. 
In addition, it includes the Aau molecular orbital. This 
model yields the largest Hamiltonian matrix among 
all the Models 1-4 (see Table |IV|. It is introduced to 
study the stability of the _R-matrix calculations with re- 
spect to extending the CAS. The scattering phase shift 
(plotted in Fig. [1]) shows encouraging agreement of the 
resonance position and width with Model 3, although the 
energy of the resonance Er — 0.169 eV is 15meV higher 
than the value obtained from Model 3. This slight shift 
towards higher energies suggests that the molecular or- 
bitals 37r„ and 4(T„ contribute more to the correlation of 
the target ground state than to the correlation of the 
(A'' -I- l)-electron system. The energy of the ground state 
of the target calculated using Model 4 is 0.66 eV lower 
than the energy calculated using Model 3 (see Table |ll|) . 
The good agreement between the Er calculated using 
Models 3 and 4 suggests that the improvement of the 
scattering results by employing the anionic molecular or- 
bitals instead of those optimized for the neutral is not 
purely an artifact due to making the representation of 
the target worse. The phase shift is rather stable with 
respect to changes of the CAS that yields different ener- 
gies of the target states. 

Since the construction of the Hamiltonian matrix Hp 
in the R-matrix calculations based on the most complex 
CAS Model 4 is computationally quite demanding (see 
Table lIVp , this was performed only for the equilibrium 
internuclear distance of the neutral target. Since those 
results show good agreement with the smaller Model 3, 
that more economical model was used to study the de- 
pendence of the ^Ilg bound and continuum states of O2 
on the internuclear distance discussed below. 

The effects of the vibrational nuclear dynamics pre- 
vent us from directly comparing the cross sections calcu- 
lated in the fixed-nuclei approximation with experimen- 
tal results for the energies of the scattering electron be- 
low 1 eV. Existence of the bound anionic state (discussed 
below) leads to well pronounced boomerang oscillations 
in the elastic scattering cross sections 0, Ulll] that do 
not appear in the fixed-nuclei calculations. Further theo- 
retical study of these effects using the energy-dependent 



vibrational frame transformation [2l|, |43| based on the re- 
sults presented here will be a subject of the forthcoming 
research. 

B. The bound and continuum state of OJ (^Tlg) 

Results of the R-matrix calculations discussed here can 
be used to study effects of the vibronic coupling in the 
electron collisions with O2 or in the photodetachment 
of O^. That research requires correct characterization 
of the bound electronic state of for a range of the 
relevant nuclear geometries in addition to the scatter- 
ing phase shifts or analytical quantum defects. Fig. [5] 
shows the analytical quantum defects for several internu- 
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FIG. 2. (Color online) Analytical quantum defects /i" (solid 
lines) calculated using Model 3 as a function of the incident 
electron energy E plotted for several internuclear separations 
Rn (a.u.). The right-hand side of Eq. [14] is plotted by the 
dashed line and the intersections with solid lines determine 
energies of the (^Hg) electronic bound state. 

clear separations calculated using Model 3 that provides 
the most physical results at i?e(02). These curves are 
smooth functions of the incident electron energy, even in 
the vicinity of the resonance. Relatively slow variation of 
f/^ with the internuclear separation _R„ makes this quan- 
tity suitable for modeling the vibrational dynamics of the 
anionic complex based on the energy-dependent vibra- 
tional frame transformation [2lj . Fig. [5] also shows the 
curve corresponding to the right-hand side of Eq. (|14p . 
Its intersection with the smooth quantum defect (if it 
exists) determines the bound-state energy of O2 (^Ilg)- 
It is well known that O2 posses one electronic 

bound state (^Ilg) with the minimum of the potential 
energy curve below the potential energy minimum of 02- 
The potential energy curves of the ground state of neutral 

02 and of the anion calculated using Model 1 and Model 

3 are plotted in Fig. [3] The electronic eigenenergies of 
the anion are calculated using Eq. (|14l) for the range of 
the internuclear separations, where the anion is bound 
and the resonance energy is taken for smaller internuclear 
distances, where the anionic state has finite lifetime. 
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FIG. 3. (Color online) Potential energy curves of the target 
O2 (^Sj) ground electronic state and ('^Ilg) state. In the 
region where the anion is not stable against autodetachment, 
the anionic curve displayed represents the real part of the cor- 
responding resonance energy. The calculations using Model 1 
and Model 3 are compared. The potential energy curves are 
plotted for each model with respect to energy of the target 
ground state at its equilibrium nuclear geometry. Model 3 
yields a fixed-nucleus electron affinity Ea ~ 0.375 eV. 



Fig.[3]shovifS that Model 1 does not yield a bound state 
of the negative ion. Although the anionic potential en- 
ergy curve crosses the the one of the neutral target at 
R » 2.51 a. u., its minimum lies above the minimum of 
the neutral potential curve. This behavior is common for 
both Models 1 and 2 and can it is explained above. The 
R-matrix study by Higgins et al. also fails to predict 
the bound state of O^. The results from the reference 
were later adjusted by shifting the i?-matrix poles to 
reproduce the experimental value of the electron affinity 
and used to study the resonant vibrational excitation of 
O2 by the electron impact Q . This sort of shift can lead 
to an inconsistency between the ff-matrix poles and am- 
plitudes that can nonphysically affect the width of the 
resonance, and has accordingly not been pursued in the 
present study. 

The i?-matrix calculations based on Model 3 that em- 
ploy the molecular orbitals optimized for the negative ion 
clearly show the bound state of and yield the elec- 
tron affinity Ea — 0.375 eV. This is in good agreement 
with the experimental value 0.448 eV [Sj] supported by 
later quantum chemical calculation by Ervin et al. [13| 
(see also Table ITTl| . The crossing point of the neutral and 
anionic potential curve obtained using Model 3 is located 
at internuclear distance 2.34 a. u., very close to the equi- 
librium geometry of the neutral target. This position is 
close to that determined in Ref. [1^ . Fig. |4] shows the 
resonance and bound-state energy of relatively to the 
energy of the neutral target ground state. The resonance 
width as function of i?„ is plotted in the inset. This vi- 
sualization allows for a direct comparison of Model 1 and 
Model 3 with Ref. [7] for multiple nuclear geometries. It 



FIG. 4. (Color online) The position Er of the ^Ilg state of 
calculated as a function of the internuclear distance R 
plotted with respect to the electronic threshold. The 

results obtained using Model 1 and Model 3 are compared 
with previously the published calculations by Higgins et al. 
[7|. The width of the resonance is compared in the inset. 

confirms that the good agreement between Model 1 and 
the reference [7] is preserved for other nuclear geometries 
than the equilibrium of the neutral target, while Model 
3 yields lower values of the resonance position and width 
as well as of the energy of the anionic bound state. 

The potential energy curves plotted in Fig. [3] and the 
analytical quantum defects plotted in Fig. [5] are the cen- 
tral physical results of this study, as well as the basis 
set configurations that seem to produce the best results. 
The resonance and bound-state energies of presented 
here show that the ab initio R-matrix calculations us- 
ing Model 3 provide O2 physically correct eigenphases 
and smooth quantum defects in the range of the collision 
energies considered here, in spite of the complicated elec- 
tronic structure, whereas the previously published ab ini- 
tio studies provide only inaccurate results. 

Another quantities essential in this context are the har- 
monic vibrational frequencies of the neutral target and of 
the bound negative ion. As one can see in Table Hill the 
values obtained using Model 1 and Model 3 are in good 
agreement with the experimental value p^ . 

VI. CONCLUSION 

The ab initio study of the electronic structure of the 
^Ilg bound and continuum state of in the approx- 
imation of the fixed nuclei is presented. The scatter- 
ing eigenphases and the analytical quantum defects are 
given as functions of the scattering energy for the range of 
the internuclear separations relevant in the resonant vi- 
brational excitation of O2 and photodetachment of O^. 
The scattering energies below 1 eV, where only one elec- 
tronic channel is open, are considered. For geometries, 
where the anionic state is not stable against electron au- 
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todetachment, the eigenphases were fitted to the Breit- 
Wigner formula (IT^ to determine the resonance posi- 
tion and width. At larger internuclear distances, where 
the anion is bound, its energy was determined from the 
smooth quantum defects using Eq. (|T4l) . All the calcula- 
tions were performed using the UK molecular ff-matrix 
program suite [H, |23| . 

The results for several different CAS models show that 
a large basis set of the CI configurations and neutral tar- 
get eigenstates is necessary in the inner region to achieve 
converged eigenphases. It is found that if the molecular 
orbitals of the neutral target are employed in the inner 
region, the convergence of the CAS model with respect 
to the number of included orbitals is too slow to provide 
physically correct characterization of the bound and res- 
onant ^Ilg state of 02^. A different method of selection 
of the CI configurations other than CAS should be used 
in that case. However, it is not straightforward to find a 
more appropriate selection scheme for the CI configura- 
tions that treats the target and the {N + iVelectron scat- 
tering complex in a balanced manner [2J, . In addi- 
tion, any such alternative selection scheme would require 
substantial changes in the existing UK _R-matrix codes. 
On the other hand, if the molecular orbitals optimized 
for the negative ion are employed, the CAS approach 
in the inner region yields the eigenphases and analytical 
quantum defects that yield the parameters of the reso- 
nance and of the bound state in good agreement with the 
experimental values. 

The problem studied here is suitable for a test of the 
alternative methods of calculation of the _R-matrix. The 
method tested here is based on the solution of the lin- 
ear system of equations individually for every scattering 
energy. This method proved more suitable than the sin- 
gle complete diagonalization, when the dimensions of Hn 
exceeds 40000, at least for typical computational hard- 
ware that is routinely available at present. It represents a 
nonperturbative alternative to the partitioned fi-matrix 
method that requires a single partial diagonalization of 
the Hamiltonian matrix Hn after which the spectrum is 
completed approximately. 
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Appendix A: Notes on the method used to solve the 
linear system 

The method of calculation of the fi-matrix from the 
Hamiltonian in the inner region Hg introduced in Se c. IIIl 
was implemented using the linear solver PARDISO [5J] 
that is based on LU factorization. This solver employs a 
sparse matrix storage scheme and it is designed to han- 



dle large and sparse matrices that cannot fully fit into 
the memory using the full storage. The UK ff-matrix 
program employs the DSYEVD subroutine from the LA- 
PACK library for the complete diagonalization. The 
DSYEVD subroutine requires storage of the full matrix 
and does not benefit from the fact that H_q is usually 
sparse. 

The sparse matrix storage is the main reason, why the 
method based on solving the system of linear equations is 
more favorable for large CAS models than the complete 
diagonalization. When the dimension of H_q is smaller 
than «40000 and when the memory is large enough to 
store it using the full storage format, the complete di- 
agonalization requires less CPU time than the solution 
of multiple linear systems for several different scattering 
energies (see Table ITV)) . Performance of both methods be- 
comes comparable above this dimension and the number 
of considered scattering energies for which the _R-matrix 
needs to be calculated decides which method requires less 
CPU time. This dimension is also close to the memory 
limit of typical contemporary computers. Comparison 
of both methods beyond this size becomes complicated 
and considerably larger CAS models can be handled only 
using methods that do not require full matrix storage. 

The CPU time required by PARDISO to solve single 
system of linear equations also depends on the density of 
Hn- In general, even very large problem can be solved 
efficiently using this subroutine, if it is sufficiently sparse. 



Appendix B: Comment on the application of the 
RMPS method 

The high number of the target states required to 
achieve a converged expansion of the total {N + 1)- 
electron wave function in the inner region naturally 
suggests that the approach employing the _R-matrix 
with pseudostates (RMPS) method might be rather ef- 
ficient. Previous studies [l^ liO, l44l . l45l | show that 
an accurate treatment of the complicated polariza- 
tion effects in the inner region requires a smaller 
number of the pseudostates than of the true eigen- 
states of the target. The RMPS method was applied 
to the problem studied here and following Gorfinkiel 
and Tennyson [3| the singly-excited configurations 
(lcrg2crg3f7glf7„2cr„l7r„)-'^*(l7rg)^(Ai)"'^ were included in 
the CI expansion of the target states, where Xi is the 
pseudocontinuum orbital Q. However, those additional 
configurations do not have any considerable contribution 
to the eigenstates of the neutral target in any of Models 
1-3, where the RMPS method was tested. As a conse- 
quence, the higher pseudostates do not improve the rep- 
resentation of the polarization effects in the scattering 
calculations. 

As an attempt to improve the representation of the 
pseudostates, the doubly excited configurations, where 
the pseudocontinuum orbital is singly occupied, were 
added. These configurations decreased the energies of 



14 



the target states. Corresponding {N + l)-electron terms 
including the pseudocontinuum orbitals were added to 
the CI expansion of the (A^ + 1) -electron wave function 
as well. However, this CI model yields an extremely large 
Hamiltonian Hn with the dimension exceeding 3 x 10^. 



The evaluation of all the matrix elements would require 
too much CPU time and the i?-matrix calculation be- 
comes intractable. Therefore, the RMPS treatment of 
the electron collisions with O2 at low incident electron 
energies should be reconsidered and it will be subject of 
future research. 
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